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Explanation Valxi^ 

albedo of soil for visible 
and near-infrared lifht 
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slope of saturation vapor pressur-^ 
curve at air temperature 
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Stef an-Boltzmann constant 5*669*10 

phase angle 
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ON THE REMOTE MEASUREMENT OF E'fAFOPATION PATES ?»OM ?APE ■VE'5’ 


SOIL UKD’^P VARIABLE CIOUD COVER 


1 . Introduction 

Periodic assessment of soil water content and evaporation rates can 
greatly benefit agriculture, hydrology, and civil works. Recent research 
has shovm that daily maximum minus minimum surface soil temperature and 
daily maximum surface soil minus air temperature can be used to estimate 
soil y/ater content and bare soil evaporation rates ( Idso et al, 1975 
a, b, c). These procedures have been developed and experimentally veri- 
fied only for clear-sky conditions, A continuous clear-sky condition is, 
however, not the normal case, particularly in the eastern US, In an ‘ 
eastern state such as liaryland, it is normal that the sky is at least 
partly cloud-covered in the afternoon. In order to investigate the 
feasibility of estimating daily evaporation rates under variable cloud 
cover, evaporation rates were calculated over 0,1-hour intervals and the 
results integrated over a 2if-hour period. The integrated value agrees 
reasonably well with calculations in which daily averages of weather 
parameters were used. The high time resolution provides a vivid history 
of the evaporation rate during the day and of the condensation rate 
during the night. 

Evaporation from soil under full or near-full cloud cover has not been 
investigated because of the lack of adequate Instrumentation for mea- 
suring thermal sky radiation and also because such cloud cover prohibits 
remote satellite measurements of reflected solar and of thermal radiation 
from soil, anyway. 


2. Experimental Procedure 


2. 1 Experimental Site 

The experiment site was at the Goddard Space Flight Center at Oreenbelt, 
Maryland, at a geographic longitude of ' ond a geographic latitude 

of 59^ 0*. T'he soil was brought there from a sugar beet field of the 

Experiment Station of the US Department of Agriculture at Peltoville, 
?-laryland, and had the following composition: 

19.6 % clay (particle sizes below 2 ) 

48.9 ■'’5 silt (particle sizes from 2 to 50 ju) 

31*5 % sand (particle sizes greater than 50 ^). 

The soil was in a hole 150 cm x 15‘9 cm large and about 25 cm deep, '^he 
site was 14 n west of an 8 n high building; thus, the soil surface was 

in the shadow of the building for about Yl^^egreeV^iTr ~ hours 
after sunrise and unshadowed thereafter until sunset. 

2.2 Tomperature Measurements 

Eight thermistors were used to obtain air temperatures at 10 and 1 50 cm 
above the soil surface and soil temperatures at 0.5, 1» ^0 and 20 cm 

below the soil surface, 

Gky radiation temperature in the band from 8 to 14 p was measured occa- 
sionally using a thermal radiometer (Fames Engineering P^T-5), "’he radio 
meter had been previously calibrated in the laboratory, using a black 
body at various known temperatures, Pecause the radiation temperature of 
a clear sky was always below the sensitivity r-^nge (-20°) of the PPT-5, 
v/e covered the aperture with .a thin foil. The foil increased the detec- 


table thermal radiation by reflecting emitted radiation from the interior 
of the sensor head back onto the detector and by emitting thermal radiation 
by itself. The configuration v/ith the foil was calibrated by viewing at 
cool objects at various temperatures (mostly the sky close to the horizon) 
and comparing the temperature readings with and without foil. 

Two pyranometers were used to measure the total radiation from the hemi- 
sphere above the horizon and the reflected radiation from the soil in the 
visible and near infrared regions of the spectrum. The up-looking pyrano- 
meter was'^placed on top of an instrument rack, 200 cm above the soil sur- 
face, whereas the down-looking pyranometer v/as suspended on a bar, 
approximately 20 cm above the soil surface. The shado’.y of the dovm-looking 
pyranometer (approximately 0,25 steradian) was usually in its own field of 
view and as a result the measured albedo v/as a few percent too low, ^ 
correction was made by dividing the pyranometer reading by the solid angle 
ratio 


T - P.Z5 
% 


j= 0,92 


2, 5 y/ind Velocity Measurements 

Two cup-aneraometers were used to obtain wind velocities at 10 and 1 50 cm 
above the soil surface, /. t wind velocities below ?. m/s, the cups v/ould 
normally not rotate, and a wind velocity of 1 rn/s v."»s assumed. 

•2.4 Soil Moisture lleesurements 

Four gypsum block soil moisture sensors were place! in the soil at depths 
of 2, 5> 10, and 20 cm. Their resistance v/as measured v/ith alternating 
current bridges. Temperature readings from thermistors at corresponding 

■ ■ --.5 


depths were used to correct gypsum block resistance for its temper-^ ture 
dependence. 

Core samples were taken from the soil and the volumetric soil moisture 
was measured gravimetric ally in order to calibrate the gypsum block rea- 
dings. 


?.» 5 Relative Humidity was measured with a sling psychrometer one to four 
times each day at 10 cm and I 50 cm above the ground. Absolute humidity f 
was comput,ed from relative humidity RH using the formula 


f = RH * f 


sat 


( 1 ) 


where is maximum v/ater content of air at temperature T (in °K) , 

'sat » «89 • ( a ) 

^e being the saturation water vapor pressure in Torr (or mm Hg) 


6755 


® = exp (-.^rt^ - 4,819 In T + 53.21). 

'X 


( 3 ) 


It was found that f (at 150 cm only) v/as usually slov/ly increasing v/ith 
time during dry periods and rapidly dropping after a rain. A linear in- 
crease of absolute humidity between consecutive readings of relative hu- 
midity v/as, therefore, assumed. Relative humidity betv/een readings was 
then computed from interpolated absolute humidity values. 


2.6 Precipitation was measured using a rain gauge. Ou-alitative notes were 
taken on percent cloud cover, haze density and v/lnd type ("gusts versus 
steady wind). 


2,7 Data Reduction 

All automatic readings v/ere taken at 6 minute intervals on a 24-channel 
recorder, readings v/ere manually entered into a Hewlett Packard model 

4 - 


9SlO desk C‘=»lcul!ator which converted the vilues to physic'll •rr'cnl tudes , 
using Inborn tory cnlibntion curves, '^he reduced d«tn were airaultnneoualy 
printed and plotted. At the seae time the cnlculntor performed n harmonic 
analysis of each curve f (t) over n T = 24-hour period, computing 43 pairs 
of Fourier coefficients, 

ax - ^ d-oo A olt 


« 


( 4 ) 


( 5 ) 


(where Id = 2 7t / T is the angular frequency and k s 0, 1,2, ,47 in- 

Lent 

7 ^ 


dicates the k-th Fourier coefficient), "’ho aj^, bj^ coefficients ware con- 


verted 


\ =7 ' 


+ b 

k 




arc t»n 


( 6 ) 
( 7 ) 




such that the original curve can be appro::ima ted by the sum of its har- 
monic components as 

The maxinum number of coefficients wns given by the number of '•vnil »bla 
storage registers in the C‘»lculntor. 


5. 


'Experimental Results 


. 1 3hor t-~.va ve solar and sky radiation 


L'ipures 1 and 2 show the shore- >ve r'>dlatior 
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iiifriired radiation on October 13/1^* 
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diys, July <?3/-9 )ctotor 15/Wi. -’r qj*.*nt dips In thy c’lr/as c'>n bo 

notic’d, dun to Incr'^'sed cloudin.'sr in th ■ •’ftcrmon. "’’ho r»'flpctod 
r.odl'itlon is nhov/n in ”i -urn J for OcloVor 13/14. 

onr-jr.yo oky rodlntton 

The thoraol rndiotlon fro'a ^ claTr-s’:y •*t">oophnre h*-3 been c-rlcul^tod 
fron the air tespentur'’ soai-onipirlol Idoo— ^nckaon 

fornul^ (Jdso and .Tackaon , 1969^ 




( 9 ) 


Assuralnj an eraisoion coefficient for wet soil of €_ ~ 0*8^, the fnctior. 

0.65 K of the thormnl sky radi-ition absorbed by the soil h'-s been c'lcu- 
a 

Intei for the October 13/14 ob.servation and Jo shown in Figure 4. Tt can 
bo seen that this energy flux is nearly con.stant. 

An attempt has been Tirade to noasuro thor.o’l sky radiation using the F'U"’-5 
r.adior.eter. Figure 5 snows the radiation temperature in the b^nd from 
8 to 14 Jim ns e function of zenith .angle on .a clear sky COctoler ?.Z ) and 
on 3 cloudy sky ( October 25 ). ''•ir temperatures on both days have been 

Q 

betw.aen 10 and 15 C. '"he r-diatlon temperature- at zenith angle from 0 
to 7ipproxlma tely 6o dsi'r • .> :s -35^C for the clo T .5ky and + 15 ^ C for the 
cloujy sky. The corre.ip niin;* emitted r ii-tion in the b tnd from 8 to 1 4 

p 

is .ani 14.3 m,V/rs. , r*sr . r ^ vely. ‘I'^houn v’ me .oured the th'^-m 1 

r!:y radiation Only in th-- '-ni from ^ to 14 |rn, th. ' c'ncluni ’n apre-TS 
inevlt-ible that the total radiation fror a cloudy .eky is sign! fic- ntly 
hirh-.'r th >n th-.t from a .ole-.r sky - nu *h ‘ the Tdsa- ^.aCKSon fO''mula fq) 
ic not -pplicablo to » clou -., r.^y. 

"h? h-.t. e. 3 ti.mito n f thfr-r.-il r i — t’an fr 'm a ci'uiiv s':v ( aurismln F = 1) 
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i^iire 3. Reflected visible and near-infrared radiation on October 13/14 
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sky it is R = 27.6 raV.'/cm . Tho difference between those two esti*mted 

SI 

fluxes, = 3.8 rn',V/cm , happens to be nenrly the sntne as the difference 
between the two measured fluxes, AH = 8.9 ra'Vcm , This result seems to 
justify our assumption that the emission coefficient of a dense cloud is 
nearly unity. 


3*3 I.ong~v/ave soil radiation 


The thermal emission from a soil at temperature T„ 

H 5 



Is 


Rg =£^<r/ ( 10 ) 

Assuming again €. = 0.35, Figure 6 shov/s the thermal emission from the 
soil as calculated for the October 15/1 4 obseryation. The energy flux is 
shown negative because it is directed away from the soil surface. 


5.4 Soil Heat Flux 


For calculating the heat flux into and out of a soil surface one usually 

has to know the thermal properties of the soil. Accordingly, an attempt 

has b^en made to compute the thermal diffusivity of the soil, 

'hi z k /p c 

from the amplitudes and ph^’ses of the temperature v/aves penetrating the 
soil as given in To.ble 1, The two appro-.ches used are described by 
Crorlslaw and Jaager (1959) 
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Figure 6. Thermal radiation of soil on October 13/14 









































on 





officlonta of those curv"3 ire listed '^•»blo I, "ho rosMlto ~ro is 
follo.v3: 


layer at 
depth 
X,- X, 

,5-1 CO 

1 -2 era 

2-5 cm 

1 

O 

n 

1 

IO->0 cm 

^in 

^from 

0.0012 

0.01 1 

0.0067 

0.0055 

0.0049 

•^2 

0.03095 

0.01 3 

0.0064 

0.oh53 

0.0059 


0.0035 

C.O 056 

0.0062 

0,0070 

0.0033 


0.0039 

0.0058 

0.0060 

0,0072 

0.0035 

aver/\se 

0.0023 

0.0059 

0.0063 

0.0065 

0.0044 


Two chnnecteristics ire noticod. First, the thorn«l dlffusivlty of the 
layer fron 0,5 to 1 cn is sirnl flc'»ntly lower thm "t layers Velow 1 cm. 
This is interpretoJ is an indication of a lower soil moisture content near 
the soil surface. Second, although the values calculated from T^ agree 
well with those cilcul'>ted from and the values calculated from ^ y 
agree well with those calculated from J^2» approaches load to 

slgnlficmtly different results, ".'hilc in the uppermost layer A/(?) equals 
only about 0, 3«(V). in the next deeper layer the situation is reversed 
and J*6(T) « 2 ^(f). 

It is believed that the diS'-ero'/ment is C'used by •’r. onerry transport 
mechanism different from heat conduc ti on, such as diffusion of W'’t*r or 
V." tO!' vapor throUf^h the soil, '>ud that the mc -dires ‘ssumlnf heat con- 
duction as the only '.raporirt rear.-r.lst ;rc not directly ?‘p''llcaile to a 
v/et soil. 


•• different arnroach .ut ti-ts 
is knOAT., f?ince Lot thermal 
••ini that of .v'-ter } - "^,16, 


fro- the '•SE'ir.ptl on th t the thermal inertia 
inertia of -n aver- c.} moist soil is p = 0 , 1 h 
ou” 'SSU'"r*.ior appf»- to re justified, '’’ho 
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hoit flux is C'»lcul-'ted using only the curve of the teuper'»turc *it the 
coil surf-co, 

Vhen the soil texnenture T -t - depth x ia n rer'.odlc function of tl*ne 
of period = 24 hours, 


T* - Tj, + & iot-mfTx ■M'i ) 


( ’3 ) 


(wherr tn = ^ is the '.TTve number of « tetjner" turo w ve whoso wevelongth 


is X ^IL , ^ bein.„ the thorm-1 diffusivity of Lhe m*»teri=»l 

'' ' Sc 


then the flux of hent Into the soil is 




( 14 ) 


** * e r s / 1 cut - 1-*^ 1/7 y i- (^^ }J 

V'e m^y use the identity 

ct B'/? Ai^ci coo ^ YT’c&o d 

sYT yOy*i^ t ^S*) 

to obtnln 

s=k^W 1 1;- rr^‘’"”!kv + 'f; +vs-<>) 

t-»l 
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/>! 
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P b’^ln.' the therml inertl'i of the soil, ••/e oV t' ■’ n 


( 15 ) 
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( 17 ) 
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‘’surface * fSvX, /uU (tof ■►•VI +*#•?•) •/•lEo'PTi 

(^ 1 ) 

s 5^ /W -f S. *) + 

( >2 ) 

.Ve aee thnt the tr^nsforintton fro"!! the Courier coefficients of the tet^- 
per'«tura curve those of flux curve follows the rules 

for the ith orsplitude ? 7 / ( 23 ) 

ond for the ith phoee • Y^, ( 24 ) 

Fifure 7 shows the soil surface temperature on October 13 / 14 . A har- 

s 

;sonic analysis results In the following expression: 

T[*C] -}ni*y «f.03 + ^ 25 ) 

Assuainp a thertr.,al inerti" typic^’l for wet soil, 

Ts 0.t>t W,*'‘ c«'* *<■' 

the heat flux into the soil surface can now be calculated, following equa- 
tions ( 23 ) and ( 24 ) 

5 ^hK^/c»h *7 102.2*)^ ( PS ) 

The first 16 pairs of fourlor coefficients of both T and S curves are 
list'.'l in Table 1 . is shov;r. in '‘ipuro 13 » 
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•’he Powen ntio 1 h nlottei in -i'ure ^U for In'? October oVnorvntlon. 

^ho ro 3 ult of o roint-py-polnt c-’lcul-'blon of " is shO'jrn in '’i^ure 15« 

The enerfy >jvoilHble for eviror-'tlon err. be errr«3s<?d by 

Th <9 evapoj'ition bscoraes neg'’tive during th<? niebt, indic'>tinif conrtonsotion 
of .vntar on zhe soil surfTCO. "ne interr-' tad flux of ©vapor-. tinf w*>ter 
during the 2ii-hour period is c^tlcul*'ted fro"i the first coefficient of eq. 


(39), 


P a • SL^OOt -/Ox /-/ X /J 


oun’.'mrv of ^'r.err'y Fluxes 


The contributions of the various kinds of energy, nveraged over 24 hours 
are sunjaerized below ( all in a.V/ca“) 


incident light 

reflected light 

ther.ml soil redi'^ticn 

theranl sky redi-tior. 

•absorbed by soil 

sent flux ttirough soil su-f'ce 

•ovnr'retion en-*rfy 

sensible neot t-"' w rr. • ij’ 
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Figure 15. Heat of evaporation (left-hand scale) and rate of evaporation (right-hand scale) 
on October 13/14. 


Goap-trison of C^lculiteJ \’v«*por-»tlon with •ubltshod Moiols 

l’sin*r as the dally averMj^a of a ••vonther par'^nater Its 0-th ^ourlor 
coefficients as plven In '^aVle i, we may calcul*’te daily evaporation 
rates from published rnodels and compare the results with our ?h-hour 
evaporation flux. 


Sink .Stren :th Approaches 
Foh?;er (1930 found empirically 


?. = 0.10 (e^ - e/) (1 ♦ 0.0071 u ) 

8 a 


( 


•Vind speed u (miles 
is shown in ^i»^uros 
speed at 1.5m, u = 


per hour) on October 15/1h et 0.1 m and 1,5 m 
16 and 17, respectively. Using the average wind 
2,4 m/c, '.ve obtain 

5.3 mrn/d. 


Four formulae given by Penman (194S) yield the follov/ing results: 

K = 0.065 - e^) u = 3.9 mm/d, ( 

K = 0.033 (6g - e^) u = 3.9 mm/d, ( 


K 


•and F 

Penman considers 


= 0.30 (1 0.014^ u) 'o - 0 .) 

(5 (1 

= 0.35 (1 * '.'> 09 '^ 1,0 '. o ^ - -• 

the last onuation to te th-* best form 


3.7 rnm/d,( 

3.1 ”’m/d. ( 
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'nrtv/y ^'••l.ancc '’■prrJacho s 


The "enmon (lOl'i) formula for ^’r. enerrv b l->nce -crro-’ch fives 
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'Isin-* a t'orrrul , of ^’.rthollc et il, for '*1 ev'^poration f’^om. 


a wot surface. 
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Figure 17. Wind speed at 150 cm on October 13/lA. 







Conclusions ^nd Recom.'nendf.tions 


The r>stim-ites of daily evapoT'^tion rates fron? nil .approaches .a^reo to 
within 5*'^ percent, periote me.asure.T.ent of evaporation rates from bare 
wet soil appears feasible, even under a lic;ht cloud cover, if no 
better accuracy is required. 

In order to improve the accuracy of ev-aporabion estimates, the follo- 
wing research is recommended; 

The evaporated mass of v/ater should be monitored quantitatively using 

4 

a weighing lysimeter, so that the accuracy of each individual .method 
can be given. 

The correlation between cloud cover and thermal sky radiation should be 
investigated quantitatively and systematically using .a radiometer which 
is sensitive in the entire range of vmvelengths from 3 to 20 jim. 

Valid general procedures can only be based on a large number of obser- 
vations at various conditions. In an approach to reduce the large num- 
ber of variables affecting evaporation rates, certain recurring v;eather 
patterns, their frequencies and characteristics should be identified, 
and one should search for possible accumulations of such patterns at 
various times in the annual cycle and at various locations in the US. 
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